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Abstract 

This  is  the  first  in  a  series  of  papers  in  which  we  present  state-of-the-art  methods  demonstrated  at  Case  for  the  estimation  of  transport  properties 
in  gas  diffusion  layers  (GDLs)  for  proton  exchange  membrane  fuel  cells  (PEMFCs).  Most  of  the  methods  used  today  for  measuring  wettability 
properties  of  GDLs  are  related  to  the  external  contact  angle  to  water.  The  external  contact  angle  however  does  not  describe  adequately  capillary 
forces  acting  on  the  water  inside  the  GDL  pores.  We  show  as  well  that  the  direct  method  of  estimation  of  the  internal  contact  angle  using  goniometry 
on  micrographs  is  impractical.  We  propose  and  describe  in  this  paper  a  method  for  estimating  the  internal  contact  angle  to  water  and  the  surface 
energy  of  hydrophobic  and  hydrophilic  gas  diffusion  media.  The  method  was  applied  to  GDLs  having  different  contents  of  hydrophobic  agent 
and  carbon  types.  The  method  can  be  applied  separately  to  different  components  of  the  GDL  including  macro-porous  substrates  and  micro-porous 
layers.  The  uncertainty  estimates  using  this  method  are  usually  within  3%  of  the  measured  value. 

©  2006  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Proton  exchange  membrane  fuel  cell;  Transport  property;  Gas  diffusion  layer;  Wettability;  Contact  angle;  Surface  energy 


1.  Introduction 

Gas  diffusion  layers  (GDLs)  are  essential  components  of 
polymer  electrolyte-based  fuel  cells,  including  hydrogen  and 
methanol  fueled  systems.  The  role  of  GDLs  within  a  proton 
exchange  membrane  fuel  cell  (PEMFC)  is  to  allow  the  gaseous 
reactants  to  move  towards  the  catalyst  layer  areas  situated  above 
the  ribs  between  the  gas  channels.  They  also  provide  a  path  for 
electrons  to  flow  between  catalyst  layers  and  bipolar  plates.  The 
GDL  also  plays  a  critical  role  in  water  management  within  the 
cell.  Cell  humidification  and  water  removal  are  both  achieved 
through  the  GDL.  GDLs  are  responsible  for  an  optimum  humid¬ 
ification  level  and  liquid  water  distribution  in  the  catalyst  layer. 
This  is  required  in  order  to  minimize  the  ohmic  resistance  of 
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the  ionomer  and  to  prevent  the  appearance  of  local  hot  spots  in 
the  catalyst  coated  membrane.  Depending  on  specific  operating 
conditions,  a  variety  of  aspects  of  the  GDL  come  into  play. 

At  the  fuel  cell  cathode,  supersaturated  conditions  are  some¬ 
times  encountered.  The  liquid  water  produced  in  the  cathode 
catalyst  layer  changes  the  effective  porosity  of  the  GDL  hin¬ 
dering  the  ability  of  oxygen  to  diffuse  towards  the  catalyst, 
a  phenomenon  referred  to  as  ‘flooding’  of  the  electrode.  This 
is  typically  the  origin  of  the  limiting  current  for  PEM  fuel 
cells.  In  addition,  resulting  occlusion  of  flow  channels  by  liquid 
water  droplets  at  the  GDL/channel  interface  have  recently  been 
implicated  in  problems  with  durability  in  cells.  Local  reagent 
starvation  seems  to  lead  to  substantial  local  voltage  excursions, 
causing  carbon  corrosion. 

The  anode  of  a  hydrogen/air  cell,  on  the  other  hand,  tends  to 
dry  out  and  it  may  be  desirable  for  the  GDL  to  act  to  retain  water. 
In  a  liquid-fed  direct  methanol  fuel  cell,  the  GDL  represents  an 
adjustable  impedance  to  reagent  delivery.  Thus,  the  wettability 
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Nomenclature 

ak  unknown  parameters  (k  =  1,2) 

A  vector  in  Eq.  (A.  10) 

C  matrix  in  Eq.  (A.  10) 

Cw  Washburn  constant  (cm5) 

/  functional  relationship  between  dependent  and 

independent  variables 
F  conditional  function 

L  weight  factor 

m  mass  absorbed  by  test  sample  (g) 

n  number  of  measurements 

R  residuals 

S  weighted  sum  of  the  squares  of  residuals 

t  time  (s) 

V  vector  in  Eq.  (A.  10) 

w  weights 

Xj  calculated  independent  variable  (J=  1,2,3) 

Xj  measured  independent  variable  (j=  1,2,3) 

y  calculated  dependent  variable 

Y  measured  dependent  variable 

Greek  letters 

y  surface  tension  (mN  m~ 1 ) 

r]  dynamic  viscosity  (mPas) 

0  contact  angle 

p  gas  density  (gem-3) 

a  standard  deviation 

Subscripts 
L  liquid 

LV  liquid-vapor  interface 

SV  solid-vapor  interface 

SL  solid-liquid  interface 

Superscripts 

d  dispersive  component  of  surface  tension 

p  polar  component  of  surface  tension 


characteristics  of  the  GDL  are  of  paramount  importance  in  its 
function  and  in  overall  fuel  cell  function. 

Given  the  importance  of  the  GDL  in  fuel  cell  function,  it  is 
somewhat  surprising  that  they  have  not  been  the  subject  of  sig¬ 
nificantly  more  work.  This  stems  partly  from  the  complexity  of 
media  themselves.  Significant  structural  complexity  underlying 
the  GDL  behavior  leads  to  difficulty  in  defining  a  set  of  easily 
measurable,  rigorously  meaningful  parameters  of  description. 

Gas  diffusion  layers  are  usually  constructed  from  macro- 
porous  substrates  coated  with  one  or  more  micro-porous  layers. 
The  morphology  of  the  GDLs  for  PEMFCs  may  vary  widely 
depending  on  the  nature  and  proportion  of  their  constituents. 
Macro-porous  substrates  for  GDLs  may  be  carbon  fiber  cloth 
(Fig.  1)  carbon  fiber  paper  (Fig.  2),  or  non-woven  carbon  mate¬ 
rials.  The  surfaces  of  the  carbon  materials  are  not  so  inert  as, 
for  example  the  pure  sp3  bonding  of  diamond  but  include  struc¬ 
tures  that  provide  some  polar  character  to  the  carbon  materials 


Fig.  1.  200  x  reflected  light  microscopy  image  of  an  ELAT®  carbon  fiber  cloth 
(E-TEK  division  De  Nora  North  America). 


used  in  GDLs.  To  increase  the  contact  angle  to  water  and  pre¬ 
vent  flooding,  macro-porous  backbones  of  gas  diffusion  layers 
are  impregnated  with  poly(tetrafluoroethylene)  (PTFE).  Micro- 
porous  layers  (amorphous  mixtures  of  carbon  and  PTFE)  are 
applied  on  macro-porous  substrates  in  order  to  improve  elec¬ 
trical  contact  with  the  catalyst  layers  and  to  add  a  further  level 
of  control  of  water  transport.  Indeed,  the  macro-porous  layer 
used  alone  will  tend  to  fill  with  liquid  water,  or  flood,  at  the  fuel 
cell  cathode.  For  a  structurally  complex  medium  such  as  this, 
defining  the  factors  controlling  wetting  is  extremely  difficult. 
However,  advancing  the  state  of  the  art  requires  development 
of  rigorous  descriptions  of  these  phenomena,  supported  by  the 
development  of  rigorous  testing  methods. 


Fig.  2.  500 x  scanning  electron  microscopy  (SEM)  image  of  a  Toray®  carbon 
fiber  paper.  The  image  shows  carbon  fibers,  clusters  of  PTFE  and  small  par¬ 
ticles  on  the  fibers  (PTFE  or  of  other  origin).  At  pore  (microscopic)  level  the 
GDL  material  may  not  be  considered  homogeneous  as  in  a  macro-homogeneous 
model.  The  internal  contact  angle  to  water  will  depend  on  the  nature  of  the  mate¬ 
rial  the  droplet  is  in  contact  with.  The  measured  internal  contact  angle  to  water 
must  represent  a  statistical  average  over  the  GDL  properties  at  microscopic  level. 
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Liquid  water  within  the  GDL  is  transported  by  capillary  and 
gas-phase  pressure  gradients  while  body  forces  such  as  gravity 
and  two-phase  drag  do  not  significantly  affect  its  momentum. 
The  properties  of  the  GDL  that  influence  the  capillary  action 
are  its  surface  energy  and  the  contact  angle  to  water.  To  support 
the  research  of  new  candidate  materials  for  GDLs  as  well  as 
for  design  considerations  one  must,  at  a  minimum,  be  able  to 
provide  quantitative  estimations  of  the  surface  tension  and  the 
internal  contact  angle  to  water  of  GDL  materials. 

Most  of  the  methods  used  today  to  estimate  the  contact  angle 
to  water  of  GDL  materials,  such  as  capillary  meniscus  height, 
goniometry  (sessile  drop)  or  Wilhelmy  (plate)  are  methods  for 
determining  the  external  contact  angle.  Usually  values  of  the 
external  contact  angle  between  120°  and  140°  or  even  higher 
have  been  reported  [1-3]  for  different  GDL  materials.  Since  the 
contact  angle  to  water  of  pure  PTFE  is  only  108°  [4]  these  large 
values  cannot  be  explained  by  the  presence  of  the  hydrophobic 
agent  inside  the  GDL  pores  but  rather  by  the  contribution  of 
GDL  surface  roughness.  During  an  external  contact  angle  mea¬ 
surement  water  does  not  experience  the  molecular  force  field 
with  the  internal  pores  of  a  porous  material,  but  the  contact  with 
a  solid  plate  having  surfaces  of  a  certain  roughness.  For  a  rough, 
structured  surface,  these  measurements  do  not  provide  a  true 
measure  of  the  interfacial  properties  of  the  material.  Instead, 
structural  aspects  dominate  the  observed  behavior  through  such 
phenomena  as  droplet  pinning.  It  is  well  known  [5]  that  surfaces 
with  a  roughness  above  approximately  1  |xm  generate  higher 
contact  angles  to  water  than  smooth  surfaces  and  a  hysteresis 
effect  when  advancing  and  receding  contact  angles  are  mea¬ 
sured.  While  in  general  the  external  contact  angle  may  describe 
only  qualitatively  the  wettability  of  GDL  materials  it  may  not 
be  used  as  a  quantitative  estimation  for  calculations  or  design 
purposes. 

In  contrast,  the  internal  contact  angle  to  water  of  gas  diffusion 
media  depends  on  the  surface  energy  at  the  pore  level.  Since  at 
the  pore  (microscopic)  level  the  structure  and  composition  of  gas 
diffusion  media  is  not  homogeneous  (see  Fig.  2),  a  meaningful 
internal  contact  angle  must  reflect  a  statistical  average  over  the 
material  properties  at  microscopic  level.  For  this  reason,  direct 
measurement  of  internal  contact  angle  using  goniometry  associ¬ 
ated  with  SEM  or  other  optical  techniques  is  impractical  given 
that  it  would  require  the  interpretation  of  a  very  large  number  of 
droplets.  Use  of  the  direct  method  for  interpretation  of  internal 
contact  angles  is  complicated  even  further  if  the  contact  point 
between  the  water  droplet  and  the  solid  phase  in  a  micrograph  is 
hidden  behind  the  droplet  or  behind  the  fiber  (see  Fig.  3).  Even 
when  this  is  not  the  case,  it  has  been  shown  [6]  that  “micro”  and 
“macro”  contact  angles  measured  using  goniometry  may  differ 
significantly. 

An  internal  contact  angle  to  water  reflecting  a  statistical  aver¬ 
age  over  diffusion  media  properties  at  the  microscopic  level  can 
be  obtained  using  the  sorption  or  Washburn  method  [7],  The 
Washburn  method  has  been  cited  in  the  past  in  the  literature 
related  to  gas  diffusion  media  [8]  as  a  possible  method  to  be 
used  to  determine  the  contact  angle  to  water,  but  no  results  or 
even  attempts  to  use  it  for  hydrophobic  GDL  materials  have  been 
reported;  this  is  not  surprising  considering  that  sorption  methods 


Fig.  3.  500  x  scanning  electron  microscopy  (SEM)  image  of  a  Toray®  carbon 
fiber  paper.  The  image  shows  water  droplets  formed  on  the  GDL  fibers.  Despite 
the  large  number  of  droplets,  all  the  boundaries  between  water  and  fibers  in  this 
image  are  situated  either  behind  the  droplets  or  behind  the  fibers,  making  the 
interpretation  of  the  contact  angle  inaccurate. 


may  not  be  used  with  hydrophobic  porous  media  when  the  test 
fluid  is  water. 

Mathias  et  al.  [3]  refer  to  the  Zisman  plot  [9]  as  a  method  to 
be  used  to  determine  the  surface  energy  of  GDL  fibers,  but  once 
again  no  results  or  attempts  to  use  it  are  reported.  Our  initial 
attempt  to  use  this  method  to  determine  the  surface  energy  of 
the  GDL  fibers  and  the  internal  contact  angle  to  water  failed. 
Since  Zisman  theory  is  a  one  component  model  for  the  surface 
energy,  it  does  not  emphasize  molecular  interactions;  for  this 
reason  it  works  best  with  non-polar  surfaces,  but  becomes  inad¬ 
equate  when  surfaces  are  even  marginally  polar  [5],  It  appears 
that  dipole-dipole  and  even  hydrogen  bonding  interactions  are 
present  between  the  gas  diffusion  media  and  the  test  fluids.  When 
such  polar  interactions  are  possible,  other  models  for  surface 
energy  must  be  employed.  It  is  beyond  the  scope  of  this  paper 
to  evaluate  the  various  surface  chemistry  models  that  have  been 
developed  and  that  are  the  subject  of  much  discussion  in  the 
literature.  See  Adamson  and  Gast  [5]  for  some  details.  From  an 
engineering  perspective,  the  two-parameter  model  of  the  sur¬ 
face  tension  provides  a  reasonable  and  practical  algorithm  for 
determining  the  internal  contact  angle.  In  this  model,  the  sur¬ 
face  tension  is  represented  by  the  sum  of  two  components:  a 
term  that  represents  the  result  of  dispersion  forces  and  a  term 
that  includes  all  other  contributions.  The  algorithm  depends  on 
curve  fitting  data  for  different  liquids  (with  different  surface 
energy  components)  sorbing  into  the  GDL  fiber  network.  In  the 
end  the  number  of  points  in  the  sample  space  of  the  fitting  is 
the  number  of  liquids  that  sorb  significantly  for  which  the  two 
components  of  the  surface  tension  have  been  estimated  by  other 
means.  Therefore,  the  number  of  degrees  of  freedom  for  the  fit 
is  a  small  number.  As  a  result  the  uncertainty  estimates  are  nec¬ 
essarily  crude  but  useful.  Note:  a  fitting  function  that  employs 
three  or  more  parameters,  for  example  the  acid-base  model  [5], 
will  not  allow  uncertainty  estimates  since  the  data  set  and  there¬ 
fore  the  degrees  of  freedom  is  too  small.  Moreover,  the  purpose 
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of  our  experiment  and  data  processing  is  to  provide  information 
that  leads  to  material  design  criteria  for  GDL  fiber  networks.  A 
discussion  of  the  fundamental  problem  of  the  molecular  level 
interpretation  of  the  many-parameter  partition  of  the  surface 
energy  or  the  equation  of  state  approach  is  deferred  for  future 
publications. 

2.  Materials  and  method 

The  technique  demonstrated  at  Case  for  estimating  the  inter¬ 
nal  contact  angle  to  water  and  the  surface  energy  of  GDL  mate¬ 
rials  uses  the  Washburn  method  [7]  for  hydrophilic  materials 
and  a  combination  of  the  Washburn  method  and  the  Owens- 
Wendt  two-parameters  theory  [4]  for  hydrophobic  materials.  In 
the  second  case,  since  liquid  water  cannot  be  used  as  test  fluid, 
the  Washburn  method  is  conducted  with  a  set  of  wetting  fluids  to 
find  their  internal  contact  angle  to  the  GDL  material.  The  Owens- 
Wendt  theory  is  used  next  to  extrapolate  the  data  obtained  with 
the  Washburn  method  and  determine  the  internal  contact  angle 
to  water  and  the  surface  energy  of  the  GDL  material.  A  least 
squares  fit  analysis  method  [  10]  is  used  to  interpret  the  data  and 
provide  error  analysis. 

Four  GDL  samples  were  prepared  by  E-TEK  division, 
DeNora  North  America  (Somerset,  New  Jersey).  The  samples 
consisted  of  a  single-sided  carbon  fiber  cloth  (the  macro-porous 
substrate  was  the  same  for  all  samples)  impregnated  with  a 
micro-porous  sub-layer  (four  combinations  of  two  different 
types  of  carbon  and  30%  and  70%  PTFE  loading,  respectively). 

The  internal  contact  angles  to  water  of  the  GDL  materials 
and  their  surface  tensions  were  determined  at  room  tempera¬ 
ture  using  a  Kriiss  Processor  Tensiometer.  In  the  tensiometer 
the  GDL  sample  is  held  by  a  metal  clamp,  which  is  attached 
to  a  sensitive  balance.  The  test  liquid  is  placed  on  a  platform 
on  a  precise  screw-type  motor,  which  raises  it  until  the  balance 
detects  the  contact  of  the  GDL  sample  to  the  liquid  surface.  A 
computer  records  the  mass  of  liquid  absorbed  by  the  GDL  as  a 
function  of  time.  The  data  acquisition  rate  was  of  10  records  per 
second. 

The  set  of  test  fluids  used  with  this  technique  and  their  sur¬ 
face  tension  components  are  summarized  in  Table  1.  All  test 
fluids  were  obtained  from  Fischer  Scientific  and  were  of  HPLC 
grade. 


Table  1 

The  test  fluids  and  their  surface  tension  with  dispersive  and  polar  components 


Test  fluids 

K £v  [mN  m“ 1  ] 

y[v  [mNm-1] 

Klv  [mN  m  1  ] 

For  Washburn  method 

w-hexane 

18.4 

0.0 

18.4 

Toluene 

26.1 

2.3 

28.4 

Acetone 

17.3 

6.4 

23.7 

Benzyl  alcohol 

30.3 

8.7 

39.0 

Diiodomethane 

50.8 

0.0 

50.8 

Methanol 

16.0 

6.7 

22.7 

For  Owens-Wendt  calculation 

Water 

19.9 

52.2 

72.1 

Source :  Kruss  processor  tensiometer  database. 


3.  Theory  and  calculation 

3.1.  The  Washburn  method  [7] 

The  Washburn  method  described  below  is  the  first  step  of  the 
analysis  and  is  used  to  determine  the  Washburn  constant  of  the 
material  and  the  internal  contact  angle  of  the  test  fluids  to  the 
GDL  material. 

The  Lucas- Washburn  Eq.  (1 )  describes  the  capillary  rise  of  a 
liquid  in  a  porous  material  when  inertia  and  gravity  forces  are 
negligible: 


v  -  ~  \  l  J 

t  CwPFKlv 

The  GDL  sample  is  first  tested  with  n-hexane  that  is  assumed 
to  have  a  zero  contact  angle  to  the  GDL  pores  (cos  0  =  1)  and 
the  Washburn  constant  of  the  material  is  determined  from  the 
Lucas- Washburn  equation.  In  Eq.  (1)  Cw  is  the  Washburn  con¬ 
stant  of  the  GDL  sample,  m  the  mass  of  liquid  absorbed  by  the 
sample  in  time  t,  /;l,  Pl  and  yix  the  liquid  viscosity,  density,  and 
surface  tension  and  0  is  the  contact  angle;  the  viscosity  and  den¬ 
sity  of  the  gas-phase  are  ignored.  The  Washburn  constant  Cw  is 
material  dependent  and  is  a  function  of  the  GDL  pore  structure. 
With  the  Washburn  constant  determined,  contact  angles  against 
at  least  four  different  wetting  liquids  (see  Table  1 )  are  measured 
using  mass-squared-versus-time  data  and  the  Lucas-Washburn 
equation  (1).  The  contact  angle  against  water  and  the  surface  ten¬ 
sion  of  the  GDL  material  are  next  calculated  from  the  previously 
determined  contact  angles  using  the  Owens-Wendt  equation. 


3.2.  The  Owens-Wendt  calculation  of  contact  angle  to 
water  and  surface  energy  [4 ] 


According  to  the  Owens-Wendt  two-parameter  model  the  sur¬ 
face  tensions  of  the  solid-vapor  and  liquid-vapor  interphases 
consist  of  two  components:  a  dispersive  one  accounting  for  van 
der  Waals  and  other  non-site-specific  interactions  and  a  polar  one 
accounting  for  dipole-dipole,  dipole-induced  dipole,  hydrogen 
bonding  and  other  site-specific  interactions  [5] 

Ksv  =  Ksv  +  Ksv  (2a) 

Klv  =  Yin  +  Yin  (2b) 

Combining  Good’s  equation  (3)  [11]  with  Young’s  equation  (4) 


Ksl  =  Tsv  +  Tlv  -  2 


TsvTlv  ' 


TsvTlv 


(3) 


KLV  cos  6  =  ysv  -  FSL 


(4) 


yields 


(1  +  COS0) 


Klv  +  Klv 


(5) 


where  0  is  the  contact  angle  between  the  test  liquid  and  the  pores 
of  the  GDL  material.  The  two  unknown  components  of  the  GDL 
surface  tension  and  y|v  in  Eq.  (5)  can  now  be  determined 
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from  the  previously  measured  contact  angles  against  the  test 
fluids  with  known  values  of  surface  tension  components 

and  Y\y  (Table  1).  A  plot  of  (1  +  cos  d)((y[y  +  Yiy)/^\[yty) 

versus  y  Klv/Klv  f°r  different  liquids  yields  the  dispersive  com¬ 
ponent  y^y  (square  of  the  y-intercept),  the  polar  component  y£v 
(square  of  the  slope)  and  consequently  the  surface  tension  of  the 
solid-vapor  interface  ysv  from  Eqs.  (2a)  and  (2b).  The  plots  for 
the  same  GDL  material  and  different  liquids  are  situated  on  a 

straight  line;  for  the  known  value  of  \J Y\y/Y\x  f°r  water,  the 
contact  angle  of  the  GDL  material  against  water  can  be  calcu¬ 
lated  from  the  ordinate  after  extrapolation. 

Alternatively,  the  adjustment  of  the  values  of  y^y  and  ygV 
in  a  least  squares  analysis  of  contact  angle  data  does  yield 
the  dispersive  and  polar  components  of  the  surface  tension  of 
the  solid-vapor  interface  and  consequently  ysv  from  Eqs.  (2a) 
and  (2b).  Moreover,  estimates  of  the  uncertainty  of  the  surface 
energy  components  and  thereby  ysv  can  be  estimated  from  the 
Hessian  matrix  computed  in  the  least  squares  procedure.  Since 
the  degrees  of  freedom  is  small,  and  all  of  the  experimental 
numbers  are  uncertain  to  a  significant  extent,  the  least  squares 
algorithm  must  be  carefully  formulated.  We  have  discussed  the 
details  in  Appendix  A. 


Table  3 

The  internal  contact  angles  against  water  and  the  surface  tensions  of  the  GDL 
materials  with  dispersive  and  polar  components 


Sample 

&h2o  1°1 

Tsv 

Tsv 

ysv 

30%  PTFE  carbon  type  1 

89  ±  3 

13  ±  1 

8  ±  2 

21  ±  2 

70%  PTFE  carbon  type  1 

101  ±  3 

13.8  ±  0.8 

3.1  ±  0.8 

17  ±  1 

30%  PTFE  carbon  type  2 

88  ±  7 

14  ±  2 

8  ±  3 

22  ±  4 

70%  PTFE  carbon  type  2 

96  ±  7 

14  ±  2 

4  ±  2 

19  ±  3 

Fig.  4.  Owens-Wendt  plot  for  the  GDL  materials  containing,  30%  PTFE  and 
carbon  type  1  and  carbon  type  2,  respectively. 


4.  Results  and  discussion 

The  average  Washburn  constants  of  the  four  GDL  materials 
and  the  internal  contact  angles  to  the  test  fluids  are  shown  in 
Table  2  (see  also  [12]).  The  Washburn  constants  were  found  to 
be  similar  among  the  four  GDL  samples  and  repeatable  with 
the  largest  standard  deviation  less  than  10%  for  the  sample  with 
30%  PTLE  and  carbon  type  2.  The  results  for  the  Washburn 
constants  are  based  on  measurements  of  five  samples  for  each 
batch.  The  results  for  the  internal  contact  angles  are  based  on 
measurements  of  three  to  five  samples  for  each  batch. 

The  internal  contact  angle  to  water  and  the  surface  tension 
of  the  four  GDL  samples  with  their  dispersive  and  polar  com¬ 
ponent  are  shown  in  Table  3.  The  margins  of  error  shown  in  the 
table  represent  80%  confidence  limit.  The  results  for  the  inter¬ 
nal  contact  angle  to  water  and  for  the  surface  tension  are  based 
on  measurements  with  four  or  five  test  fluids  for  each  batch. 
Typically  uncertainties  within  3°  are  obtained  for  the  internal 
contact  angle  when  five  test  fluids  are  used.  Larger  uncertain¬ 
ties  were  obtained  for  the  samples  containing  carbon  type  2  for 
which  only  four  test  fluids  were  used.  The  results  show  that 


the  materials  with  the  same  PTLE  content  have  similar  internal 
contact  angles  to  water  and  free  surface  energies.  As  expected, 
the  GDL  materials  containing  70%  PTLE  have  a  higher  internal 
contact  angle  to  water  and  lower  surface  tensions  than  the  mate¬ 
rials  containing  30%  PTLE.  All  measured  contact  angles  have 
values  less  than  the  contact  angle  of  pure  PTLE  to  water.  All 
samples  have  similar  values  for  the  dispersive  component  of  the 
surface  tension.  The  polar  component  of  the  surface  tension  is 
similar  for  the  samples  with  same  PTLE  content;  it  is  the  polar 
component  which  differentiates  between  the  samples  contain¬ 
ing  30%  and  70%  PTLE,  respectively.  This  is  the  reason  why 
the  initial  attempts  to  use  the  one  parameter  Zisman  model  for 
surface  tension  failed. 

Ligs.  4  and  5  present  the  Owens-Wendt  plots  for  the  samples 
obtained  with  four  and  five  testing  fluids. 

Remark.  While  the  results  shown  here  were  obtained  for 
GDL  materials  containing  both  macro-porous  substrate  and 
micro-porous  layer,  we  have  demonstrated  this  method  at  Case 
for  macro-porous  substrates  separately,  for  micro-porous  layers 
applied  on  solid  substrates  and  for  catalyst  layers. 


Table  2 

The  GDL  Washburn  constants  and  the  internal  contact  angles  against  the  test  fluids 

Sample  Washburn  constant  (x  10-5  cm5)  Internal  contact  angle,  6  [°] 


Toluene 

(CH3)2CO 

Benzyl  alcohol 

ch2i2 

MeOH 

30%  PTFE  carbon  type  1 

1.45  ±  0.03 

47  ±  3 

21  ±  2 

68  ±1 

89.9 

0 

70%  PTFE  carbon  type  1 

1.33  ±  0.06 

55  ±  1 

48  ±  3 

71  ±  5 

89.4 

47  ±2 

30%  PTFE  carbon  type  2 

1.4  ±  0.1 

43  ±  3 

31  ±  4 

0a 

89.7 

0 

70%  PTFE  carbon  type  2 

1.47  ±  0.04 

50  ±  1 

39  ±  4 

84±4a 

89.8 

42  ±4 

a  Value  not  used  in  subsequent  analysis. 
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Fig.  5.  Owens- Wendt  plot  for  the  GDL  materials  containing  70%  PTFE  and 
carbon  type  1  and  carbon  type  2,  respectively. 


5.  Conclusions 


Appendix  A 

An  alternative  to  the  graphical  Owens- Wendt  method  is  the 
method  of  least  squares  [10],  which  provides  in  addition  the 
standard  deviation  of  the  estimated  parameters.  While  there  are 
computer  programs  that  will  do  a  least  squares  analysis  of  data 
fit  to  a  model,  the  number  of  degrees  of  freedom  of  the  GDL 
contact  angle  is  sufficiently  small  that  the  weights  must  be  care¬ 
fully  estimated  and  uncertainties  included  in  details  as  we  outline 
next.  The  fundamental  assumption  is  that  the  measured  quan¬ 
tities  obey  a  Gaussian  distribution;  the  law  of  large  numbers 
generally  assures  that  this  is  the  case,  see  for  example  Mann 
[10]  and  citations  therein.  However,  clearly,  the  number  of  liq¬ 
uids  that  are  suitable  is  limited  and  therefore  the  degrees  of 
freedom  in  determining  the  fits  is  limited.  It  is  therefore  nec¬ 
essary  to  explain  in  some  detail  the  algorithm  that  we  used  to 
provides  both  the  estimate  of  parameters  and  an  estimate  of  the 
uncertainty  in  their  values. 

The  physical  law  (5)  is  rewritten  as 


The  external  contact  angle  to  water  is  not  indicative  of  the 
capillary  forces  acting  on  the  water  inside  the  pores  of  GDL 
materials.  The  external  contact  angle  provides  information  on 
the  wetting  phenomena  on  the  external  surfaces  of  the  GDL, 
while  the  internal  contact  angle  provides  information  on  the 
wetting  phenomena  in  the  GDL  pores.  The  former  depends  on 
the  surface  roughness  more  than  on  the  material  composition 
(the  measured  contact  angle  is  outside  the  range  of  the  contact 
angles  measured  for  each  individual  component),  while  the  latter 
depends  on  the  pore  structure  and  materials  at  the  pore  level. 

It  is  not  practical  to  measure  the  internal  contact  angle  to 
water  using  goniometry  on  micrographs.  Since  at  the  micro¬ 
scopic  level  the  GDL  materials  are  not  homogeneous,  a  method 
for  the  estimation  of  internal  contact  angles  must  provide  contact 
angles  representing  a  statistical  average  over  the  GDL  properties 
at  the  pore  level. 

We  have  demonstrated  and  presented  a  method  for  estimation 
of  the  internal  contact  angle  to  water  and  of  the  surface  tension 
of  hydrophobic  and  hydrophilic  GDL  materials.  The  method 
combines  the  Washburn  technique  with  the  Owens- Wendt  two- 
parameter  theory  and  the  method  of  least  squares.  The  method 
was  applied  to  four  GDL  materials  containing  two  different  car¬ 
bon  types  and  30%  and  70%  PTFE,  respectively.  The  samples 
containing  70%  PTFE  have  higher  contact  angles  to  water  and 
lower  surface  tensions.  All  measured  contact  angles  to  water  are 
less  than  the  contact  angle  to  water  of  pure  PTFE.  Uncertain¬ 
ties  within  3°  degrees  for  the  contact  angle  are  usually  obtained 
with  this  method  when  at  least  five  test  fluids  are  employed. 
The  method  is  applicable  to  macro-porous  substrates  and  micro- 
porous  layers  of  GDL  materials  as  well  as  to  catalyst  layers. 
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y  —  +  /a  2*2  (A.  1 ) 

or  as  a  functional  relationship  between  the  dependent  and  inde¬ 
pendent  variables 

y  =  /(at,  a2,  x\,  X'i)  (A.2a) 

where  the  dependent  variable  is 

y  =  (Klv  +  KlvXcos6>+  !)  (A  2b) 

the  independent  variables  are  x\  —  y£ v,  X2  —  y£v  and  the 
unknown  parameters  are  a\  —  y[?v  and  a 2  =  yfv. 

The  residuals  are  defined  as  the  difference  between  the  mea¬ 
sured  and  calculated  values  of  the  dependent  and  independent 
variables,  respectively 


Ryt  =  Yi  —  yt,  RXji  =  Xjj  -  xjj,  7  =  1,2  (A.3) 

The  weights  are  defined  as  the  reciprocal  of  the  squares  of  the 
measurement  uncertainties  (standard  deviations) 


(A.4) 


Notice  that  ayi  must  be  computed  from  Eq.  (A. 2b),  since  we  can 
estimate  the  uncertainties  of  [y^y,  y£ y,  $}  from  experimental 
information.  A  simple  calculation  with  Eq.  (A.2b),  Mann  [10] 
provides  the  following  formula  for  the  error  propagation 


2 


o 


(A.4b) 
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this  research. 


and  thereby  wyi  can  be  estimated.  Also  recognize  that  since  the 
uncertainty  in  the  { yf'y .  yj\, }  must  also  be  included  through  vuyi . 

The  method  of  least  squares  consists  in  determining  the  val¬ 
ues  of  the  parameters  a\  and  a 2  which  minimize  the  weighted 
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sum  of  the  squares  of  residuals 


S  =  E 


"■vX,  + 

;=  i 


(A. 5) 


This  amounts  to  a  maximum  likelihood  algorithm  for  a  Gaussian 
distribution  function  of  noise  involved  in  measuring  the  various 
parameters. 

The  conditional  function  is  defined  as 


F‘  =  yi  -  f(ai,a2,  x\,x2)  (A.6) 

The  initial  guess  of  the  unknown  parameters  is  a®  and  a®-  The 
initial  value  of  the  conditional  function  is  estimated  from  mea¬ 
surements  and  the  initial  guess  of  the  unknown  parameters 

F{,  =  Y,  -  f(a(l  4  X\ ,  X2)  (A. 7) 


The  weight  factor  for  each  measurement  is  defined  as 


The  optimized  values  of  the  unknown  parameters  are  calculated 
as 


A  =  C~lV 


(A. 9) 


where 


A  = 


a®  —  a\ 

a®  —  a2 


C  = 


A  1  dF ‘  dF1  A  1  dF‘  3 Fl 
2—f  Lj  da\  da i  /—'Ll  3fli  da2 

i=l  i=l 

^  1  3 Fl  dF'  A  1  dF'  dF' 
2—f  Lj  da\  da2  /—'Ll  da2  da2 

i=l  i=l 


V  = 


J/  1  3  F‘  , 

^Ti~dlnFo 


i=i 


^  1  dF'  , 

^Ti~dJ,Fo 


L  i=l 


(A.  10) 


The  standard  deviations  for  the  unknown  parameters  are  esti¬ 
mated  by 


J“k 


n  —  2 


1/2 


(q*V/2 


(A.  11) 


The  classical  algorithm  [10]  for  determining  a\  and  a2  is  the 
following 


1.  Collect  the  data  measurements  F;,  X\ ,  X2  and  the  standard 
deviations  cryi ,  aXj . ; 

2.  Make  the  initial  guess  for  the  unknown  parameters  a®  and 

3.  Using  the  initial  guess  for  the  unknown  parameters  and 
the  measurements,  calculate  the  conditional  function  (Eq. 
(A. 7))  and  its  partial  derivatives  from  (Eq.  (A.6)); 

4.  Calculate  the  weight  factors  (Eq.  (A. 8)); 

5.  Calculate  the  elements  of  matrix  C  and  vector  V  (Eq. 
(A.  10)); 

6.  Invert  matrix  C; 

7.  Calculate  vector  A  (Eq.  (A.  10))  and  a\  and  a2\ 

8.  Use  newly  calculated  a \  and  a2  as  the  initial  guesses  for 
step  2; 

9.  Continue  steps  2-8  until  the  change  in  a\  and  a2  is  less  than 
a  prescribed  value; 

10.  Calculate  the  standard  deviations  for  a\  and  a2  (Eq. 
(A.  11)); 

References 

[1]  C.  Lim,  C.Y.  Wang,  Electrochim.  Acta  49  (2004)  4149. 

[2]  D.  Wood,  J.  Davey,  F.  Garzon,  P.  Atanassov,  R.  Borup,  Proceedings  of 
the  206th  Meeting  of  the  Electrochemical  Society,  Hawaii,  Honolulu, 
2004. 

[3]  M.  Mathias,  J.  Roth,  J.  Fleming,  W.  Lehnert,  in:  W.  Vielstich,  H.A.  Lamm, 
A.  Gasteiger  (Eds.),  Handbook  of  Fuel  Cells  Fundamentals,  Technology 
and  Applications,  vol.  3,  Wiley,  2003,  pp.  1—21. 

[4]  D.K.  Owens,  R.C.  Wendt,  J.  Appl.  Polym.  Sci.  13  (1969)  1741— 
1747. 

[5]  A.W.  Adamson,  A.P.  Gast,  Physical  Chemistry  of  Surfaces,  6th  ed.,  Wiley, 
New  York,  1997. 

[6]  H.  Giesche,  in:  E.  Schiith,  K.S.W.  Sing,  J.  Weitkamp  (Eds.),  Handbook 
of  Porous  Solids,  vol.  1,  WILEY-VCH  Verlag  GmbH,  2002,  pp.  309- 
349. 

[7]  E.W.  Washburn,  Phys.  Rev.  18  (1921)  273-283. 

[8]  C.Y.  Wang,  in:  W.  Vielstich,  A.  Lamm,  H.A.  Gasteiger  (Eds.),  Handbook 
of  Fuel  Cells  Fundamentals,  Technology  and  Applications,  vol.  3,  Wiley, 
2003,  pp.  337-347. 

[9]  W.A.  Zisman,  ACS  Advances  in  Chemistry  Series  #43  (1961), 
1-51. 

[10]  J.A.  Mann  Jr.,  in:  E.  Matijevic,  R.J.  Good  (Eds.),  Surface  and  Col¬ 
loid  Science,  vol.  13,  Plenum  Press,  New  York,  1984,  213  pp.; 
Also  see  J.R.  Wolberg,  Data  Analysis  using  the  Method  of  Least 
Squares:  Extracting  the  Most  Information  from  Experiments,  Springer, 
2006. 

[11]  R.J.  Good,  L.A.  Girifalco,  J.  Phys.  Chem.  64  (1960)  561. 

[12]  M.J.  Bluemle,  M.Sc.  Thesis,  Case  Western  Reserve  University,  Cleveland, 
OH,  2004. 


